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Abstract 

Background:  Sodal  stress  mouse  model,  based  on  the  resident  intruder  paradigm  was  used  to  simulate  features 
of  human  post  traumatic  stress  disorder  (PTSD).  The  model  involved  exposure  of  an  intruder  (subject)  mouse  to  a 
resident  aggressor  mouse  followed  by  exposure  to  trauma  reminders  with  rest  periods.  C57BL/6  mice  exposed  to 
SJL  aggressor  mice  exhibited  behaviors  suggested  as  PTSD  in  mouse  phenotypes:  intermittent  freezing,  reduced 
locomotion,  avoidance  of  the  aggressor  associated  cue  and  apparent  startled  jumping.  Brain  tissues  (amygdala, 
hippocampus,  medial  prefrontal  cortex,  septal  region,  corpus  striatum  and  ventral  striatum)  from  subject  (aggressor 
exposed:  Agg  E)  and  control  C57Biy6  mice  were  colleaed  at  one,  10  and  42  days  post  aggressor  exposure  sessions. 
Transcripts  in  these  brain  regions  were  assayed  using  Agilent's  mouse  genome  wide  arrays. 

Results:  Pathways  and  biological  processes  associated  with  differentially  regulated  genes  were  mainly  those 
thought  to  be  involved  in  fear  related  behavioral  responses  and  neuronal  signaling.  Expression  based  assessments 
of  activation  patterns  showed  inaeased  activations  of  pathways  related  to  anxiety  disorders  (hyperactivity  and  fear 
responses),  impaired  cognition,  mood  disorders,  circadian  rhythm  disruption,  and  impaired  territorial  and  aggressive 
behaviors.  In  amygdala,  activations  of  these  pathways  were  more  pronounced  at  earlier  time  points,  with  some 
attenuation  after  longer  rest  periods.  In  hippocampus  and  medial  prefrontal  cortex,  activation  patterns  were 
observed  at  later  time  points.  Signaling  pathways  associated  with  PTSD  comorbid  conditions,  such  as  diabetes, 
metabolic  disorder,  inflammation  and  cardiac  infarction,  were  also  significantly  enriched.  In  contrast,  signaling 
processes  related  to  neurogenesis  and  synaptic  plasticity  were  inhibited. 

Conclusions:  Our  data  suggests  activations  of  behavioral  responses  associated  with  anxiety  disorders  as  well  as 
inhibition  of  neuronal  signaling  pathways  important  for  neurogenesis,  cognition  and  extinction  of  fear  memory. 
These  pathways  along  with  comorbid  related  signaling  pathways  indicate  the  pervasive  and  multisystem  effeas  of 
aggressor  exposure  in  mice,  potentially  mirroring  the  pathologic  conditions  of  PTSD  patients. 
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Background 

Post-traumatic  stress  disorder  (PTSD)  is  an  anxiety  dis¬ 
order,  which  can  result  from  exposure  to  a  traumatic  event 
(s).  PTSD  can  manifest  itself  in  a  number  of  ways  such  as  a 
persistent  re-experiencing  of  the  traumatic  event,  arousal 
avoidance,  numbing,  fear,  and/or  starde  responses  [1]. 
PTSD  patients  may  also  exhibit  impaired  cognitioa  work¬ 
ing  memory  deficit,  mood  disorders  and  diminished  social 
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activities  [2].  Common  PTSD  co-morbid  conditions  include 
depression,  obesity,  diabetes,  peripheral  inflammation,  car¬ 
diovascular  problems,  and  metabolic  disorders  [3]. 

Anxiety  disorders,  including  PTSD,  encompass  hetero¬ 
geneous  conditions  that  are  associated  with  abnormalities 
in  the  fear  response  circuitry  in  brain  structures  such  as 
the  medial  prefrontal  cortex,  amygdala,  hippocampus,  and 
nucleus  accumbens  (Additional  file  1:  Figure  SI).  Changes 
in  the  activity  and  functional  connectivity  of  these  brain 
regions  lead  to  abnormalities  in  the  perception  and  inter¬ 
pretation  of  traumatic  events  and  to  the  development  and 
maintenance  of  anxiety  disorders  [4-7].  For  example,  the 
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hippocampus  (HC)  is  implicated  in  contextual  fear  learn¬ 
ing,  in  trauma  memory  consolidation,  and  retrieval  (intru¬ 
sive  re-experiencing  of  the  traumatic  event);  the  amygdala 
(AY)  is  associated  with  cue  conditioning,  hyper-vigilance, 
heightened  arousal,  learning,  and  expressing  fear  behavior 
[8,9];  the  medial  prefrontal  cortex  (MPFC)  is  involved  in 
cognition,  emotional  regulation,  and  fear  extinction;  and 
the  nucleus  accumbens  (NAc,  part  of  the  ventral  striatum, 
VS)  is  associated  with  reward-related  behavioral  abnor¬ 
malities  such  as  depression  [10].  Other  brain  regions,  such 
as  corpus  striatum  (ST)  and  septal  region  (SE),  are  impli¬ 
cated  in  processing  and  responding  to  reward  and  aversive 
stimuli,  lack  of  motivation,  and  recall  of  intrusive  trau¬ 
matic  cues  in  the  form  of  flashbacks  [11,12]. 

Impaired  extinction  of  fear-potentiated  startle  and  en¬ 
hanced  cue  conditioning  in  these  brain  regions  (of  trau¬ 
matized  patients  and  animal  models)  can  be  manifested 
in  the  forms  of  hypervigilance,  arousal,  and  avoidance 
symptoms  [13].  In  particular,  failure  to  inhibit  fear  learn¬ 
ing  and  fear  memory  consolidation  are  considered  to  be 
precipitating  factors  in  the  development  of  PTSD  [14]. 
With  regards  to  cognition,  PTSD  patients  showed  signifi¬ 
cantly  less  activation  in  sensory  association  areas,  suggest¬ 
ing  diversion  of  attention  from  the  presented  stimuli, 
perhaps  due  to  increased  focus  on  the  elicited  trauma 
memory  [15-17].  It  has  been  suggested  that  cognitive  im¬ 
pairments  exhibited  by  people  with  PTSD  result  from  in¬ 
trusive  flashback  memories  that  transiently  interfere  with 
ongoing  cognitive  processing  [17]. 

Molecular  and  cellular  encoding  of  traumatic  events  and 
behavioral  responses  are  presumed  to  be  reflected  in  synap¬ 
tic  plasticity  or  changes  in  the  activity  and  functional  con¬ 
nectivity  among  AY,  HC  and  mPFC  [18-21].  Specifically, 
changes  in  signaling  molecules  associated  with  synaptic 
transmission  and  plasticity  are  implicated  as  a  primary  sub¬ 
strate  for  fear  learning  and  memory,  thus  putatively  leading 
to  PTSD  in  human  and  to  behavioral  features  of  PTSD  in 
animals.  This  could  occur  because  the  expression  levels 
and  types  of  signaling  molecules  determine  spike-timing- 
dependent  plasticity  in  relation  to  traumatic  perceptions 
and  responses.  This  could  lead  to  either  a  long-term  synap¬ 
tic  potentiation  (LTP)  increase  in  synaptic  strength  and  in¬ 
crease  in  excitatory  post-synaptic  potential,  potentially 
enhancing  the  manifestation  of  traumatic  disorder,  or  to 
long-term  synaptic  depression  (LTD)  leading  to  a  decrease 
in  synaptic  strength  and  a  decrease  in  EPSP  size,  suppres- 
sesing  fear  extinction  processes  [22,23].  Understanding  the 
molecular  underpinnings  of  fear  learning  and  memory  fol¬ 
lowing  trauma  exposure  or  presentation  of  an  aversive 
stimulus  would  be  critical  to  identifying  prognostic  bio¬ 
markers  and  for  developing  an  intervention  to  apply  post¬ 
trauma,  before  the  onset  of  pathological  symptoms. 

Subsequent  to  trauma  exposure,  the  symptoms  of 
stress  reaction  typically  develop  over  varying  amounts  of 


time  (typically  from  days  to  months),  resulting  in  hetero¬ 
geneous  pathologies.  The  heterogeneity  of  PTSD  symp¬ 
toms  suggests  that  its  etiology  is  diverse,  and  there 
are  still  no  known  or  accepted  molecular  biomarkers 
for  diagnosis  of  PTSD.  Current  pharmacotherapies  for 
PTSD  are  applicable  after  symptoms  manifest,  and  pri¬ 
marily  consist  of  selective  serotonin  reuptake  inhibitors. 
There  are  no  FDA-approved  pharmacological  interven¬ 
tions  available  for  the  treatment  of  traumatized  individ¬ 
uals  to  forestall  the  development  of  PTSD  prior  to  the 
onset  of  pathological  symptoms.  Many  factors  hinder 
advances  towards  identifying  prognostic  biomarkers  and 
potential  pre-symptomatic  therapeutic  targets  for  PTSD. 
One  of  these  factors,  despite  great  advances  in  brain 
imaging,  is  a  lack  of  real  time  visualization  of  cellular 
reactions  in  specific  brains  regions  of  individuals  at 
risk  for  PTSD  and  patients  with  the  disorder.  Hence 
animal  models  are  required  to  obtain  appropriate  brain 
regions  for  genomic,  genetic  and  other  high  throughput 
molecular  studies. 

Increasing  evidence  suggests  that  precipitating  factors 
such  as  chronic  stress  (including  social  defeat  stress)  in¬ 
duce  changes  in  the  functional  connectivity  within  the 
fear  circuitry,  and  such  changes  mediate  trauma-related 
behavior  alerations  in  animal  models  of  anxiety  disorders 
[24,25].  Such  behavioral  responses  include  startle  re¬ 
sponse,  hyperactivity,  avoidance,  freezing,  grooming, 
rearing,  and  deficiencies  in  novel  object  recognition  and 
in  territorial  behavior.  The  molecular  and  cellular  origins 
of  these  behavioral  abnormalities  have  been  suggested  to 
include  changes  in  glutamatergic  and  GABAergic  synap¬ 
tic  plasticity,  dopamine  neuron  excitability,  epigenetic 
and  transcriptional  mechanisms,  and  neurotrophic  factors 
[26-28].  However,  the  genome-wide  molecular  basis  for  the 
interplay  of  different  behaviors  in  anxiety  disorders  is  not 
very  clear.  We  hypothesize  that  chronic  exposure  to  stress 
alters  gene  expression  patterns  which  are  important  for  the 
functional  activity,  connectivity  and  signaling  among  neu¬ 
rons  of  these  brain  regions,  thus  leading  to  the  develop¬ 
ment  of  anxiety  and  depression-like  behaviors. 

We  and  others  have  shown  that  aggressor  exposed 
(Agg-E)  social-stress  mouse  models  elicit  behaviors  such 
as  submissive  posture,  freezing  during  locomotion,  vertical 
rearing,  grooming,  aggressor  barrier  avoidance,  hyper¬ 
activity,  jumping,  and  impaired  cognition  (deficiency  in 
novel  object  recognition)  [29,30].  These  might  be  consid¬ 
ered  to  be  PTSD-like  phenotypes. 

In  this  study,  our  objective  was  to  correlate  our  previ¬ 
ous  behavioral  observations  [29]  with  changes  at  the 
molecular  level  in  brain  regions  critical  for  fear  learning, 
fear  memory  consolidation/extinction  and  fear  response/ 
expression.  Here  we  attempted  to  understand  how  envir¬ 
onmental  stimuli  (exposure  to  aggressor)  and  transcrip- 
tome  interact  and  influence  each  other  in  the  context  of 
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behaviors  suggesting  fear.  Such  approaches  are  integral 
to  our  understanding  and  treatment  of  stress -related  dis¬ 
orders  such  as  PTSD  (interaction  of  genes,  brain,  and 
behavior).  Although  a  candidate  gene  approach,  in  con¬ 
junction  with  endophenotyping,  would  be  important  to 
solve  the  puzzle  of  PTSD  etiology,  genome-wide  screen¬ 
ing  for  transcriptome  changes  has  the  potential  to  lead 
to  new  candidate  targets  and  pathways  specific  to  stress 
disorders,  and  to  facilitate  or  hasten  the  one  gene  or  few 
candidate  genes  approaches.  It  may  also  reveal  new  mo¬ 
lecular  mechanisms  and  connections,  which  may  not  be 
apparent  from  a  single  gene  or  genotyping  approach. 

Toward  this  goal,  we  performed  global  gene  expression 
profiling  of  different  brain  regions  implicated  in  fear  and 
anxiety  processing.  Genome- wide  transcriptome  changes 
were  assessed  in  brain  regions  collected  from  male  C57BL/ 
6  mice  exposed  to  male  SJL  aggressor  mice.  The  behaviors 
elicited  suggest  a  possible  PTSD-in-mouse  phenotype: 
immobility,  avoidance  of  aggressor-associated  stimuli 
(aggressor  odor),  jumping  (indicating  hyperactivity),  freez¬ 
ing  during  locomotion  (suggesting  fear),  and  reduced  loco¬ 
motion.  Transcriptome  alterations  in  six  brain  regions: 
HC,  AY,  MPFC,  VS  (contains  NAc),  SE,  and  ST  (Additional 
file  1:  Figure  SI)  from  Agg-E  mice  were  measured  using 
Agilent  s  mouse  genome-wide  arrays  to  identify  transcripts 
and  associated  pathways  that  reveal  potential  molecular 
mechanisms  of  stress  disorders.  Profiles  of  genome-wide 
transcriptome  changes  were  carried  out  at  different  time 
points  to  assess  how  the  time  course  dynamics  of  tran¬ 
scripts  indicate  molecular  events  associated  with  traumatic 
fear  learning  and  memory  along  the  developmental  trajec¬ 
tories  of  traumatic-anxiety  disorders. 

We  identified  a  number  of  differentially  expressed  genes 
(DEGs)  across  different  brain  regions  and  time  points. 
Functional  and  pathway  analyses  of  DEGs  suggested  pos¬ 
sible  roles  in  anxiety-related  behavioral  responses,  synaptic 
plasticity,  neurogenesis,  inflammation,  obesity  (metabolic 
disorder),  and  cardiac  infarction.  In  particular,  transcripts 
important  for  the  synthesis  of  neurotransmitters,  gener¬ 
ation  and  development  of  the  hippocampus,  axonogenesis, 
dendritic  branching,  splice  variant  processing,  and  dopa¬ 
minergic  and  serotonergic  pathways  were  affected. 

Results 

Behavioral  evaluation  of  mouse  model 

We  established  [29]  a  rodent  model  manifesting  PTSD- 
like  behavioral  features.  We  believe  that,  because  the  stres¬ 
sor  depends  on  antagonistic  social  interaction  within  a 
conspecific,  this  model  offers  potential  advantages  over 
the  brief  foot-shock  model  in  simulating  traumatic  events 
of  conflict-like  situations.  Especially  as  both  the  stressor 
and  the  behaviors  assessed  are  essentially  within  the  social 
realm.  Concerns  [31]  that  the  repeated  social  stress  model 
may  allow  the  animals  to  habituate  rather  than  triggering 


the  PTSD-like  features  was  countered  by  implementing 
the  ‘randomness'  in  the  occurrence  of  the  life-threatening 
conditions.  Thereby  we  incorporated  a  critical  dimension 
of  traumatic  psychosocial  stress  with  regard  to  the  “etho- 
logical  validity”  of  the  unpredictable  and  uncontrollable 
nature  of  the  PTSD  simulating  trauma.  The  hypothesis 
was  further  validated  as  we  found  the  mice  displaying 
many  traumatic  features  after  a  long  time  interval.  The 
‘habituation'  related  concern  was  further  mitigated  as  the 
stressed  mice  displayed  significantly  reduced  urine  mark¬ 
ings,  a  signature  of  inhibited  territorial  behavior  in  the 
course  of  10-day  stress  session.  Territorial  urine  markings 
of  Agg-E  mice  were  three-fold  lower  than  those  of  control 
mice  with  linear  regression  significance  of  p  =  0.02  [29] . 
Post  10-day  Agg-E,  stressed  mice  maintained  decreased 
but  not  significant  urine  marking  (p  =  0.07)  compared  to 
the  controls.  During  the  same  time  course,  we  observed 
an  increased  bodyweights  of  Agg-E  mice.  Linear  regres¬ 
sion  model  showed  significant  difference  (p<  0.001)  in 
gain  of  body  weights  by  the  stressed  cohorts.  In  addition, 
Agg-E  groups  showed  significantly  increased  temperature 
(p  =  0.006),  and  30%  decreased  but  not  statistically  signifi¬ 
cant  corticosterone  levels  ((p  =  0.06). 

In  addition,  the  present  model  supports  a  set  of  prede¬ 
fined  validity  criteria  for  an  animal  model  of  PTSD 
[31,32].  For  instance,  the  “face  validity''  of  symptoms  simu¬ 
lating  PTSD-like  features  was  manifested  by  introducing 
the  stressed  mice  to  the  contextual  cues  immediately  after 
withdrawing  the  trauma  as  well  as  at  a  delayed  interval  or 
rest  period  (the  latter  being  equivalent  to  nearly  three 
years  of  human  life).  We  observed  a  significant  display  in 
associative  fear  memory,  anxiety  and  hypo-  and  hyper¬ 
responses  immediately  after  the  trauma-withdrawal;  and 
many  of  these  behavioral  traits,  such  as  the  restrained  tail 
rattling  (an  agonistic  behavior,  hypo-response)  and  groom¬ 
ing  (an  anxiety-like  feature,  hyper-response)  sustained 
through  the  delayed  interval.  Specifically,  while  the  num¬ 
ber  of  control  animals  displaying  tail  rattling  gradually 
increased  with  the  passage  of  42  days,  none  of  the  stressed 
mice  displayed  tail  rattling  even  42  days  after  the  end 
of  stress  exposure.  Moreover,  the  duration  of  grooming 
displayed  by  the  stressed  mice  remained  2.5  fold  higher 
(p  <  0.05)  than  the  controls,  after  42  days  post- 10  days  of 
Agg-E.  Likewise,  freezing,  a  hypo-active  fear  response 
remained  significant  from  day  1  (p  <  0.0001)  to  28  days 
(p  <  0.05)  post- 10  day  stress.  Suppressing  a  natural  instinct 
namely  the  vigilance  of  novel  object,  the  stressed  animals, 
showed  significant  avoidance  of  the  aggressor's  partition 
(p  <  0.0001)  1  day  after  stress;  however  Agg-E  mice  exhib¬ 
ited  decreased  avoidance  at  42  days  of  rest  (delayed) 
period  post  10  day  exposure  session. 

The  present  model  attested  an  essential  criterion  of 
PTSD  namely  the  persistence  (or  re-emergence)  of  a  con¬ 
textual  fear  memory.  We  further  noticed  a  ‘dose-response' 
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in  displaying  PTSD  feature,  as  in  comparing  the  S-day  and 
lO'day  long  traumatic  burdens,  the  later  proved  more 
deleterious.  The  “construct  validity”  of  the  rodent  PTSD 
model  representing  the  cellular  and  molecular  processes 
was  underlined  by  the  atypical  increase  of  many  blood 
cells'  load,  previously  associated  with  stress  [33].  Signifi¬ 
cant  myocardial  atrophies  caused  by  this  model  persisted 
after  the  delayed  interval,  which  potentially  implied  the 
cardiac  damages  caused  by  PTSD  primarily  on  the  war 
veterans.  Additional  histological  analysis  showed  lack  of 
the  expected  increase  in  dendritic  spine  density  of  pyram¬ 
idal  neurons  in  MPFC  (medial  prefrontal  cortex,  an  essen¬ 
tial  brain  region  for  fear  memory  extinction  and  memory 
sorting).  Together  this  model  strategically  modified  the  so¬ 
cial  stress  model  to  reasonably  simulate  aspects  of 
combat-related  PTSD,  and  validated  some  of  the  essential 
criteria  of  PTSD  in  rodent-models. 

Gene  expression  analyses  of  brain  regions 

Feature-extracted  and  quantile  normalized  microarray 
data  from  brain  regions  were  analyzed  to  assess  stress - 
affected  biological  processes  and  pathways  in  Agg-E 
mice  as  compared  to  controls. 

Differentially  expressed  genes  (DEGs),  compared  between 
Agg-E  and  control  (C-ctrl)  mice  for  each  brain  region  (AY, 
HC,  MPFC,  SE,  ST  and  VS),  were  identified  at  four  differ¬ 
ent  time  points  (T5R1,  T5R10,  TlORl  and  T10R42) 
(Table  1)  using  Moderated  T-test,  p  <  0.05  (Figure  1).  For 
the  10-day  Agg-E  groups,  there  were  more  DEGs  after 
the  longer  rest  period  (T10R42)  except  for  ST  which 
showed  a  large  decrease  in  DEGs  as  compared  to  the 
TlORl  group.  The  numbers  of  DEGs  of  5-day  Agg-E 
groups  were  greater  at  T5R1  for  AY,  HC  and  SE  regions; 
at  T5R10,  the  numbers  of  DEGs  increased  for  MPFC, 
ST  and  VS,  while  by  contrast  the  numbers  were  greatly 
decreased  for  HC  and  SE. 

The  greater  divergence  between  10-day  Agg-E  and  C-ctrl 
groups  after  42  days  of  rest  may  be  due  to  either  of  two 
factors.  One  factor  may  be  due  to  the  faster  recovery  of 


the  C-ctrl  mice  from  handling,  confinement  and  hunger 
stressors  compared  to  the  rate  of  recovery  of  the  Agg-E 
groups.  Another  factor  may  be  that  multiple  processes  lead¬ 
ing  to  a  return  to  homeostasis  versus  further  consolidation 
of  PTSD-like  disorders  (or  combination  of  both)  may  in¬ 
volve  expression  of  many  more  and  different  sets  of  genes 
in  the  Agg-E  mice. 

Functional  and  pathway  enrichments  and  visualization  of 
associated  networks 

Functional  and  pathway  enrichments  of  DEGs  of  the  dif¬ 
ferent  brain  regions  led  to  the  identification  of  modular- 
networks  (functional  interaction  networks)  consisting  of 
genes  implicated  in  anxiety- related  behavioral  responses 
and  the  underlying  synaptic  processes  (Figure  2).  Each 
module  of  this  network  has  many  more  member  DEGs, 
which  were  excluded  for  clarity  (visualization)  to  show 
only  representative  members.  Also,  the  list  of  pathway 
and  functional  modules  identified  were  longer  than 
shown  in  Figure  2  (Additional  file  2:  Table  SI).  Member 
DEGs  involved  in  long-term  memory,  associative  learn¬ 
ing,  and  limbic  system  development  were  largely  from 
the  HC;  DEGs  associated  with  startle  and  fear  responses 
were  largely  from  the  AY;  and  DEGs  associated  with  cir¬ 
cadian  rhythm,  cognition,  neurogenesis,  locomotory  be¬ 
havior,  dopaminergic  and  serotonergic  pathways  were 
from  different  brain  regions. 

Separation  of  each  module,  adding  more  member 
DEGs,  and  coloring  nodes  based  on  their  expression 
level  in  HC  and  AY  for  the  one-day  post  10-day  Agg-E 
group  (TlORl),  showed  how  a  particular  behavioral  re¬ 
sponse  or  biological  process  behaves  in  these  two  im¬ 
portant  brain  regions  (Figures  3  and  4).  For  example, 
modular  network  DEGs  involved  in  associative  learning 
were  largely  induced  in  AY,  whereas  DEGs  associated 
with  aggressive  and  territorial  behaviors  were  largely 
suppressed  in  AY,  suggesting  involvement  in  impaired 
territorial  and  aggressive  behavior  (Figure  3).  This  infer¬ 
ence  is  also  corroborated  by  the  finding  of  significant 


Table  1  Summary  of  experimental  conditions  and  collected  tissues 


Study  group 
(time  points) 

Experimental  and  control 
groups  (5  mice/group) 

^Subject  mouse  strain 

Aggressor  strain 

Exposure 
length  (days) 

Post  Agg  E 
length  (days) 

T5R1 

Agg  E  T5R1 

C57BL/6 

SJL 

5 

1 

C  ctr  5R1 

C57BL/6 

SJL 

5 

1 

T5R10 

Agg  E  T5R10 

C57BL/6 

SJL 

5 

10 

Cctr  T5R10 

C57BL/6 

SJL 

5 

10 

TlORl 

Agg  E  TlORl 

C57BL/6 

SJL 

10 

1 

Cctr  TlORl 

C57BL/6 

SJL 

10 

1 

T10R42 

Agg  E  T10R42 

C57BL/6 

SJL 

10 

42 

C  ctr  T10R42 

C57BL/6 

SJL 

10 

42 

^collected  brain  regions:  HC:  hippocampus,  AY:  amygdala,  MPFC:  medial  prefrontal  cortex,  VS:  ventral  striatum  or  nucleus  accumbens,  SE:  septal  region, 
ST:  corpus  striatum. 
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Figure  1  Differentially  expressed  genes  (DEGs)  across  brain  regions  at  different  time  points  (T5R1,  T5R10,  T10R1  and  T10R42).  The 

scatter  plots  show  both  log2  fold  changes  and  negative  logio  p  values  in  six  brain  regions  (across)  and  four  time  points  (down);  the  numbers  of 
up  and  down  regulated  transcripts  in  each  group  are  also  shown.  Key:  amygdala  (AY),  hippocampus  (HC),  medial  prefrontal  cortex  (MPFC),  septal 
region  (SE),  corpus  striatum  (ST),  and  ventral  striatum  (VS),  T:  number  of  days  of  trauma  or  aggressor  exposure;  R:  post  trauma  tissue 
collection  days;  Agg  E;  aggressor  exposed. 


reduction  of  urine  marking.  Many  of  the  DEGs  associ¬ 
ated  with  circadian  rhythm  were  suppressed  in  AY,  sug¬ 
gesting  a  circadian  rhythm  disruption,  whereas  those 
involved  in  corticotrophin  releasing  hormone  signaling 
were  largely  induced  in  AY  (Figure  4),  indicating  aggres¬ 
sion  (trauma)  potentiated  hypothalamic-pituitary-adrenal 
(HPA)-axis  signaling,  thus  leading  to  anxiety-related  be¬ 
havioral  responses.  The  expression  patterns  of  some  of 
the  important  nodes  showed  different  directions  in  AY 
and  HC.  The  differential  responses  between  these  two 
important  brain  regions  are  consistent  with  suboptimal 
communication  between  the  putative  fear  response  center 
(AY)  and  the  center  associated  with  contextual  processing 
(HC).  Examining  fewer  member  genes  among  those 


associated  with  some  of  stress-induced  processes  and 
pathways  such  as  glucocorticoid  receptor  signaling, 
neurotransmitter  secretion,  inflammation  and  growth 
factor  receptors,  we  still  see  a  mix  of  same  and  opposing 
directions  of  expression  of  these  genes  between  AY  and 
HC  (Additional  file  1:  Figure  S2).  Some  of  the  DEGs  in¬ 
volved  in  social  withdrawal,  immobility,  long-term  memory, 
and  startle,  anxiety  and  fear  responses  were  down-regulated 
in  AY  (Additional  file  1:  Figure  S3).  Many  of  these  down- 
regulated  genes  were  also  found  to  be  associated  with 
neurotrophic  factors  signaling,  fear  extinction,  and  func¬ 
tions  related  to  emotional  regulation. 

As  suppressed  or  induced  expression  of  associated 
genes  does  not  necessarily  correspond  to  the  activation 
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Figure  2  Modular  networks  of  DEGs  associated  with  circadian  clock,  behavioral  response  and  synaptic  transmission;  each  module 
(nodes  of  the  same  color)  forms  a  functional  module. 


or  inhibition  state  of  the  corresponding  behavioral  re¬ 
sponses,  we  calculated  an  activation  z-score  for  each  sig¬ 
nificantly  enriched  behavior  and  biological  process  based 
on  the  expression  values  of  associated  DEGs  in  each 
brain  region  at  each  time  point.  This  approach  is  im¬ 
portant  to  further  identify  which  behavioral  responses 
and  biological  processes  may  have  been  activated  or 
inhibited.  For  example,  Figure  5  shows  the  activation 
pattern  of  anxiety-like  behavior  in  AY,  HC  and  MPFC  at 
both  early  (T5R1,  TlORl)  and  late  (T5R10,  T10R42)  rest 
periods.  Here,  it  is  more  obvious  how  the  expression  dir¬ 
ection  of  each  DEG  contributes  to  the  predicted  activa¬ 
tion  (red  central  node)  or  inhibition  (blue  central  node) 
of  the  anxiety-like  behavior  (Figure  5).  The  pattern  be¬ 
came  clearer  not  only  within  a  single  network  but  also 
in  comparison  among  different  brain  regions  and  time 
points.  At  TlORl,  genes  associated  with  anxiety- like  be¬ 
havior  are  activated  in  AY  but  inhibited  in  the  MPFC, 
which  is  consistent  with  reported  observations  of  exagger¬ 
ated  activation  of  AY  and  delayed  regulation  by  MPFC 
during  early  stage  amygdala-dependent  fear  responses 
[34,35].  On  the  other  hand,  at  longer  rest  periods  (T5R10, 
T10R42),  the  opposite  pattern  of  activation  occurred  be¬ 
tween  AY  and  MPFC  (Figure  5),  which  may  indicate 
that  the  attenuation  of  the  anxiety-like  response  in  AY 
was  probably  due  to  emotional  regulation  by  MPFC. 
Anxiety-like  behavior  also  showed  activation  or  attenuation 


at  earlier  and  at  later  rest  periods,  respectively,  in  AY 
(Additional  file  1:  Figure  S4).  Activation  z-scores  for 
fear  conditioning  and  startle  response  in  AY,  HC  and 
MPFC  also  showed  different  activation  or  attenuation 
patterns  across  time  points  and  regions  (Additional  file  1: 
Figures  S5  and  S6). 

The  activation  z-score  for  each  enriched  process  at  every 
time  point  showed  increased  activation  of  many  of  the 
anxiety-related  behavioral  responses  at  one-day  post  Agg-E 
(particularly  in  AY  at  TlORl)  and  a  slight  attenuation  (but 
still  an  activated  pattern)  at  longer  rest  periods  (Figure  6). 
A  similar  approach  for  visualizing  the  signaling  pathways 
related  to  synaptic  plasticity  and  neurogenesis  showed 
largely  inhibited  patterns  (Figure  7). 

Whereas  we  observed  increased  predicted  activation  of 
pathways  implicated  in  behavioral  responses,  there  was, 
in  contrast,  a  decreased  predicted  activation  in  a  broad 
spectrum  of  pathways  pertaining  to  synaptic  plasticity 
and  neurogenesis  (Additional  file  1:  Figure  S7).  There 
were,  however,  a  few  important  exceptions  such  as  loss 
of  neurons,  long-term  synaptic  depression  (activated  in 
AY,  HC  Sc  MPFC)  at  T5R1,  and  activation  of  many  of 
these  signaling  functions  in  AY  at  TlORl,  and  in  HC 
(and  to  some  extent  in  MPFC)  at  T10R42  (Additional 
file  1:  Figure  S7).  This  may  be  related  to  many  of  the  ele¬ 
vated  neuronal  activities  of  AY  since  it  is  the  center  of 
fear  and  anxiety  responses  at  one-day  post  Agg-E.  On 
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AY  T10R1  HC  T10R1 


Figure  3  Networks  of  DEGs  involved  in  associative  learning,  fear  memory  extinction,  and  aggressive  and  territorial  behaviors. 
Nodes  are  colored  using  AY  TlORl  and  HC  T10R1. 


the  other  hand,  after  the  longer  rest  period  (T10R42), 
neuronal  signaling  may  be  more  active  in  the  HC  due  to 
long-term  memory  consolidation  and  related  processes 
and  in  the  MPFC  since  it  is  involved  in  regulating  emo¬ 
tional  responses. 

Our  assays  and  analyses  showed  that  many  genes  re¬ 
lated  to  PTSD-associated  behaviors  such  as  startle  re¬ 
sponse  (Figure  6  and  Additional  file  1:  Figure  S7),  rearing, 
cognitive  impairment,  hyperactivity,  anxiety,  addiction, 
avoidance  and  fear  conditioning  (Additional  file  1: 
Figure  S7),  had  differential  activation  across  tissues  and 
time  points,  the  larger  activation  generally  being  at  one- 
day  post- 10  days  of  Agg-E  (TlORl)  sessions  (Figure  6 
and  Additional  file  1:  Figure  S7).  Synaptic  transmissions, 
including  secretion  and  concentration  of  neurotransmit¬ 
ters,  long-term  potentiation,  and  synaptogenesis  (Figure  7) 
were  largely  inhibited  at  one-day  post- 10  days  of  Agg-E 
(TlORl)  and  10-days  post-5  days  of  Agg-E  (T5R10). 
We  observed  that  neurotransmission  and  synaptic  plasti¬ 
city  tend  to  be  inhibited  across  tissues  and  time  points 


compared  to  behavioral  responses  and  neurological  disor¬ 
ders,  indicating  impaired  cognition  and  mood  disorders 
(Figure  6),  which  are  important  facets  of  anxiety  disorders. 

Networks  of  transcripts  important  in  behavioral  stress 
responses,  synaptic  transmission  and  circadian  rhythm 
Modular  networks  of  DEGs  from  different  brain  regions 
have  been  significantly  associated  with  associative  learn¬ 
ing,  long-term  memory,  and  regulation  of  serotonergic 
and  dopaminergic  synaptic  transmissions,  cognition,  neuro¬ 
genesis,  limbic  system  development,  startle  and  fear  re¬ 
sponses,  and  circadian  rhythm  (Figure  2).  These  pathways 
and  biological  processes  are  important  molecular  signaling 
implicated  in  traumatic  responses  leading  to  anxiety-related 
disorders. 

Comparative  transcriptome  changes  among  brain  regions 
important  in  stress  response  circuitry 
Transcriptome  alterations  in  AY  and  HC  were  compared 
for  overlapping  and  differential  biological  processes  and 
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Figure  4  Networks  of  DEGs  involved  in  circadian  rhythm  and  corticotropin  releasing  factor  signaling.  Nodes  are  colored  using  AY  T10R1 
and  HC  TlORl. 


pathways.  These  include  a  number  of  transcripts  re¬ 
ported  to  be  involved  in  startle  and  fear  responses,  social 
withdrawal,  long-term  memory,  circadian  rhythm,  cor¬ 
ticotropin  releasing  hormone  receptor  signaling,  associa¬ 
tive  learning,  aggressive  and  territorial  behavior,  and  fear 
memory  extinction  (Figures  3  and  4  and  Additional  file  1: 
Figure  S3).  Red  nodes  show  activation  and  blue  nodes  show 
inhibition  of  these  particular  behaviors.  Many  showed  op¬ 
posite  regulation  in  AY  compared  to  their  expression  in 
HC  for  the  one-day  post  10-days  Agg-E  session  (TlORl) 
(Figures  3  and  4  and  Additional  file  1:  Figure  S3).  These 
data  suggest  that  there  are  opposing  forces  acting  in 
AY  compared  to  those  molecular  events  in  HC  under 
severe  stress  or  traumatic  conditions  possibly  with  incipi¬ 
ent  stress  disorders. 

Common  and  specific  transcripts  for  5-  and  10-day 
aggressor  exposure  sessions  with  different  post  Agg-E 
tissue  collection  rest  (delay)  days 

Among  DEGs  from  brain  regions  at  different  time 
points,  2,470  were  associated  with  behavioral  stress  re¬ 
sponses,  synaptic  plasticity  and  transmissions,  circadian 
rhythm,  obesity  and  diabetes,  and  inflammatory  re¬ 
sponse.  Of  the  2,470  DEGs,  62  were  specific  to  TlORl, 
106  to  T5R1,  103  to  T5R10,  401  to  T10R42,  and  306 


were  common  among  all  time  points.  The  306  common 
DEGs  were  significandy  (Bonferroni,  <7  <0.05)  associated 
with  axon  guidance,  neurotransmitter  secretion,  opioid 
signaling,  G-protein  mediated  events,  downstream  TCR- 
signaling,  regulation  of  synaptic  plasticity,  associative 
learning,  synapse  assembly,  CNS  neuronal  axonogenesis, 
long-term  potentiation,  long-term  depression,  activation 
of  NMDA  receptor  and  synaptic  vesicles  exocytosis 
(Additional  file  1:  Figure  S8).  The  T5Rl-specific  106 
DEGs  were  associated  with  positive  regulation  of  chronic 
inflammatory  response,  associative  learning,  dendritic 
spine  development,  negative  regulation  of  axon  exten¬ 
sion,  and  the  regulation  of  production  of  interleukins  10 
and  12.  Thel03  DEGs  specific  to  T5R10  were  associated 
with  regulation  of  the  perception  of  pain,  response  to 
carbohydrate  stimuli,  and  organic  acid  transport.  The 
TlORl -specific  62  DEGs  were  associated  with  positive 
regulation  of  circadian  sleep/wake  cycle,  regulation  of 
nerve  growth  faictor  receptor  signaling  pathways,  nega¬ 
tive  regulation  of  cardiac  muscle  cell  apoptosis,  and 
cytotoxic  activity;  and  the  DEGs  specific  to  T10R42 
were  associated  with  inflammatory  and  adaptive  immune 
responses.  It  appeared  that  some  of  the  biological  pro¬ 
cesses  and  pathways  relevant  to  traumatic  stress  re¬ 
sponses  started  occurring  early  (during  T5R1).  On  the 
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other  hand,  most  of  the  DEGs  found  only  after  the  lon¬ 
ger  rest  periods  were  associated  with  pain  perception 
(T5R10)  and  the  immune  response  (T10R42). 

Activation  patterns  of  behavioral  response  and  comorbid 
conditions 

Traumatic  stress  disorder-related  behaviors  such  as  fear 
conditioning,  anxiety-like  behavior  and  startle  responses 
showed  differential  activation  patterns  across  the  brain 
regions  of  AY,  HC  and  MPFC  at  different  time  points 
(Figure  6  and  Additional  file  1:  Figure  S7).  Fear  conditioning- 
related  genes  were  activated  in  the  HC  and  MPFC  (except 
T5R10)  throughout  the  experimental  time  courses,  but  in 
the  AY,  this  was  activated  at  the  early  defeat  period  and 
attenuated  at  longer  defeat  periods  (Figure  6).  Early  induc¬ 
tion  of  transcripts  implicated  in  fear  conditioning  sug¬ 
gested  activation  of  this  behavioral  response  early, 
consistent  with  the  observation  that  AY  responds  first  to 
traumatic  stress.  In  the  case  of  anxiety-like  behaviors, 
transcripts  associated  with  this  process  were  induced  in 
MPFC  after  the  longer  rest  periods  (T5R10  and  T10R42), 
which  is  an  indication  of  a  delay  in  MPFC  responses  to 


anxiety-like  disorders  compared  to  AY  (Figure  5).  On  the 
other  hand,  transcripts  associated  with  anxiety  show 
that  AY  responds  at  early  periods  and  this  response  was 
greatly  diminished  with  longer  rest  period  after  the 
traumatic  event  (Additional  file  1:  Figure  S4).  Although 
signs  of  inhibition  of  anxiety-like  behavior  at  longer  rest 
periods  in  AY  indicated  a  sense  of  recovery,  its  activa¬ 
tion  in  HC  and  MPFC  in  particular  (Figure  5)  suggested 
that  it  may  precipitate  in  the  forms  of  mood  disorder 
and  impaired  cognition. 

Comorbid  conditions  of  anxiety  disorders  such  as  in¬ 
flammation,  obesity,  diabetes  and  cardiac  infarction  were 
also  among  significantly  enriched  processes  and  path¬ 
ways  (Table  2).  These  observations  suggest  the  pervasive 
and  debilitating  effect  of  resident-intruder  social  stress 
on  the  health  of  the  Agg-E  mice. 

Transcript  -  brain-region  interactome  across  time  points 
We  constructed  transcript  -  brain-region  interactome 
using  transcripts  that  passed  p-value  less  than  0.01  and 
more  than  2  fold  across  brain  regions  and  time  points 
(Figure  8).  Persistent  with  what  we  observed  with  whole 
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Table  2  Differentially  regulated  genes  (DEGs)  significantly  associated  with  different  behavioral  responses,  neuronal 
signaling  processes  and  comorbid  conditions 

q  Value  Associated  DEGs 


Behavioral  response 

anxiety 

4.37E  17 

anxiety  like  behavior 

3.10E  10 

cognition 

9.54E  24 

cognitive  impairment 

2.94E  10 

conditioning 

4.06E  12 

depressive  disorder 

7.59E  10 

emotional  behavior 

3.83E  15 

learning 

1.25E  21 

locomotion 

6.56E  1 1 

long  term  memory 

1.37E  15 

memory 

5.63E  20 

Mood  Disorders 

2.80E  10 

neurological  signs 

2.63E  1 1 

neuromuscular  disease 

1.45E  10 

object  recognition  memory 

1.69E  09 

post  traumatic  stress  disorder 

5.32E  14 

recognition  memory 

2.36E  13 

spatial  learning 

3.72E  10 

spatial  memory 

2.95E  1 1 

Neuronal  signaling  and  neurogenesis 

activation  of  neurons 

1.22E  10 

atrophy  of  neurons 

4.63E  10 

axonogenesis 

1.84E  11 

cell  death  of  hippocampal  neurons 

1.51  E  09 

concentration  of  corticosterone 

1.28E  13 

concentration  of  cyclic  AMP 

3.58E  09 

degeneration  of  neurons 

2.07E  09 

dendritic  growth/branching 

6.77E  1 1 

depolarization  of  cells 

1.34E  10 

disorder  of  basal  ganglia 

2.28E  10 

excitation  of  neurons 

4.23E  12 

length  of  dendrites 

1.78E  12 

long  term  potentiation 

4.96E  15 

long  term  potentiation  of  synapse 

1.44E  12 

loss  of  dorsal  root  ganglion  cells 

4.01  E  10 

mobilization  of  Ca2+ 

3.95E  10 

neuronal  cell  death 

1.68E  15 

APP,BDNF,CCK,CNR1,CREB1,CRH,CRHR1,DRD2,  FKBP5,FOS,MAPK1,NPY,NTRK2 
APP,BDNF,CCK,CRH,CRHR1  ,MAPK1 

A^l ,APP,BDNF,CNR1  ,CREB1  ,CRH,CRHR1 ,  DRD2,EP300,ERK1/2,FOS,GNAI  1 , 
HDAC3,KF,MAPK1  ,MAPK3,NGF,NPY,NR3C1  ,NTRK1  ,NTRK2,PI3K  (complex),THRB 

APP,BDNF,CNR1,CRH,DRD2,FEP,NGF,NPY 

APP,BDNF,CNR1  ,CREB1  ,CRHR1  ,DRD2,ERK1  /2,  GNA1 1  ,HDAC2,MAPK1  ,NPY 

APP,BDNF,CCKBR,CREB1  ,CRH,DRD2,FKBP5,  NR3C1  ,NTRK1  ,THRB 

APP,BDNF,CCK,CNR1  ,CREB1  ,CRH,CRHR1  ,DRD2,  MAPK1  ,MAPK3,NGF,NPY,NTRK2 

A^l ,APP,BDNF,CNR1  ,CREB1  ,CRH,CRHR1  ,DRD2,  EP300,ERK1/2,FOS,GNAI1  ,HDAC3, 
HOMER1  ,MAPK1  ,MAPK3,NGF,NPY,NR3C1  ,NTRK1  ,NTRK2,  PI3K  (complex), THRB 

APP,BDNF,CCKBR,CNR1  ,CRH,CRHR1  ,DRD2,FEP,  NGF,NR3C1  ,STAT3 

APP,BDNF,CREB1,CRH,DRD2,ERK1/2,GNAI1,  HDAC3,MAPK1,NTRK2 

A^l ,APP,BDNF,CNR1  ,CREB1  ,CRH,CRHR1  ,DRD2,  EP300,ERK1/2,GNAI1  ,HDAC3, 
H0MER1,MAPK1,NGF,NPY,NR3C1,NTRK1,NTRK2,  PI3K  (complex) 

Aiai  ,APP,BDNF,CCKBR,CREB1  ,CRH,DRD2,  FKBP5,NR3C1  ,NTRK1  ,PFCG1  ,THRB 

A^l ,APP,BDNF,CCKBR,CNR1  ,CREB1 ,  CTNNBl ,  DRD2,FOS,FEP,NGF,NPY,NTRK2,PRF 

Aiai ,APP,BDNF,CCKBR,CNR1  ,CREB1 ,  CTNNBl ,DRD2,F0S,FEP,NGF,NPY,NR3C1 , 
NTRK2,  PRF 

APP,CNR1,EP300,HDAC3,NGF 

ADRAl  A,ADRA1  B,ADRA2A,ADRA2B,ADRA2C,  ADRB1  ,CNR1  ,CNR2,DRD1  ,DRD2, 
DRD3,DRD4,DRD5,HTR1A,NR3C1,SFC6A4 

APP,CNR1,EP300,GNAI1,HDAC3,MAPK1,NGF,  NR3C1 

APP,BDNF,CREB1  ,CRHR1  ,DRD2,FOS,NGF,NTRK2 

A^1,APP,BDNF,CNR1,CREB1,DRD2,NGF,NPY,  NR3C1,PI3K  (complex) 

APP, BDNF,CCK,CNR1, CRH, FOS,FEP,NGF,PI3K  (complex) 

APP,BDNF,NGF,NTRK1  ,NTRK2 

A^l  ,APP,BDNF,CCK,CREB1  ,DRD2,FOS,NGF,NTRK1  ,NTRK2 
Aiai  ,APP,BDNF,CRHR1  ,DRD2,FEP,NGF,  NTRK2 
APP,CNR1  ,CRH,CRHR1  ,DRD2,FEP,NPY,NR3C1 ,  NTRK2,STAT3 
ADCY,APP,BDNF,CCK,CNR1,CRH,DRD2,FEP,NGF,  PRF 
Aiai ,APP,BDNF,CNR1  ,CREB1  ,EP300,NGF,NTRK2,  STAT3 
Aiai ,APP,BDNF,CNR1  ,CRH,CRHR1  ,DRD2,NGF,  NTRK2 
APP,BDNF,CCK,CRH,FOS,FEP,NPY 

Aiai,APP,BDNF,CCKBR,CNR1,CREBl,  CTNNBl, DRD2,FOS,FEP,NGF,NPY,NTRK2,PRF 

APP,BDNF,CCK,CRH,FOS,NGF,NPY,NTRK2 

APP,BDNF,CRH,CRHR1  ,NGF,NTRK2 

APP,BDNF,CNR1  ,CREB1  ,CRHR1  ,DRD2,ERK1  /2,  GNAI 1  ,HDAC2,MAPK1  ,MAPK3,NGF, 
NPY,NTRK2,  PI3K  (complex) 

APP,BDNF,CNR1,CREB1,DRD2,GNAI1,HDAC2,  MAPK3,NTRK2 
BDNF,NGF,NTRK1,NTRK2 

BDNF,CCK,DRD2,GNAI1,GNB2,  inositol  triphosphate,NGF,NPY,PI3K  (complex), 
PFCG1,RXFP3,SRC 

Aiai  ,APP,BDNF,CNR1  ,CREB1  ,CRHR1 , CTNNBl  ,DRD2,EP300,ERK1/2,FOS,HDAC3, 
FEP,MAPKl,NGF,NR3Cl,NTRKl,NTRK2,PI3K(complex),SRC,STAT3 
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Table  2  Differentially  regulated  genes  (DEGs)  significantly  associated  with  different  behavioral  responses,  neuronal 
signaling  processes  and  comorbid  conditions  (Continued) 


neuroprotection 

4.24E  09 

A^l ,APP,BDNE,CNR1  ,CRH,CRHR1  ,LEP,NGE,PI3K  (complex) 

neurotransmission 

3.60E  1 7 

APP,BDNE,CCK,CCKBR,CNR1  ,CRH,CRHR1 ,  CTNNBl  ,DRD2,ERK1  /2,GN Al  1  ,HDAC2, 
LEP,NGE,  NPY,NR3C1,SRC 

proliferation  of  neuronal  cells 

2.13E09 

BDNF.CNRI, CTNNBl, DRD2,ERK1/2,LEP,NGF,  NTRK2,PI3K  (complex), STATS 

quantity  of  ACTH  in  blood 

4.32E  09 

CNR1,DRD2,NR3C1,NTRK2 

quantity  of  neurons 

4.71  E  12 

APP,BDNE,CREB1,DRD2,ERK1/2,LEP,NGE,NPY,  NTRK1,NTRK2,PRL,THRB 

quantity  of  pituitary  cells 

2.10E  11 

CRH, CTNNBl, DRD2,LEP,PRL,THRB 

release  of  acetylcholine 

3.84E  09 

APP,BDNE,CNR1  ,CRH,DRD2,NGE 

release  of  Ca2+ 

3.82E  09 

AI^1,APP,BDNE,CCK,DRD2,  inositol  triphosphate,LEP,NGE,PI3K  (complex), 
PLCG1,SRC 

release  of  catecholamine 

3.32E  14 

APP,BDNE,CNR1,CRHR1,DRD2,LEP,NGE,NPY,  NTRK1,NTRK2,PRL,SRC 

release  of  dopamine 

1.84E  11 

APP,BDNE,CNR1,CRHR1,DRD2,LEP,NGE,NTRK2,  PRL 

release  of  L  glutamic  acid 

1.12E  11 

ADCY,APP,BDNE,CCK,CNR1,CRH,DRD2,LEP,NGE,  NPY 

secretion  of  catecholamine 

2.06E  12 

A^l  ,APP,CNR1  ,CRH,CRHR1  ,DRD2,NPY,NTRK1 

secretion  of  corticosterone 

4.27E  12 

BDNE,CRH,CRHR1,DRD2,LEP,NPY 

secretion  of  neurotransmitter 

3.44E  10 

CCK,CNR1  ,CRH,DRD2,NGE,NPY,NTRK1  ,NTRK2 

stimulation  of  neurons 

4.28E  13 

APP,BDNE,CCK,CRH,EOS,LEP,NGE,NPY,NTRK2 

synaptogenesis 

2.22E  10 

APP,BDNE,CREB1 , CTNNBl ,HDAC2,NGE,NTRK1 ,  NTRK2 

synthesis  of  neurotransmitter 

7.69E  10 

Aiai  ,APP,BDNE,CRH,NGE,NTRK1 

synthesis  of  steroid 

3.71  E  09 

APP,BDNE,CREB1  ,CRH,CRHR1  ,LEP,MAPK1  ,NR3C1  ,PI3K  (complex),PRL,PROXl 

Comorbid  conditions 

chronic  inflammation 

1.56E  11 

ACE,ADORA2B,ADORA3,AGER,AOC3,C3,CCL5,CCR5,CD24A,CD28,CD47,CNRl, 
C0X2,EDNRA,EABP4,ECER1  A,ECER1  G,ECGR1  ,ECGR3,GIE,HSPD1  ,ID01  ,IL1  B,ITGA2, 
JAK2,LBP,LTA,MIE,PDE5A,PIK3CG,PLA2G4A,PRKCA,PTGER3,PTGS2,STAT5A,STAT5B, 
TAC1,TGM2,TLR4,TNERSE11A,TNERSE1A,TNESE11,H2  DMA,H2  DMB1,H2BEM, 
HADH,HNE1  A,HPSE,HSPD1  ,H^,ICA1  ,IENGR2,IKBKG,IL1 2B,IL1  R1  ,IL2,ILDR2,INPP5K, 
INPPL1,INSR,IRAK1,IRE1,IRS3,JAK2,KCNJ11,KHK,KRAS,LARS,LEP,LEPR,LEPR0T,LIPE, 
LMAN1L,LPIN1,LTA,  MAP3K7,MAPK8,MYD88,MY05A,NCK1,NEKB1,NEKB2,NEKBIB, 
NEKB1,NEKB2,NEKBIB,RELA 

myocardial  infarction 

6.49E  12 

Aiai ,Aiai ,APP,CNR1  ,CREB1 , CTNNBl ,DRD2,EP300,  ERK1/2,E0S,HDAC2,HDAC3, 
KL,LEP,MAPK1,  MAPK3,NPY,NR3C1,STAT3,THRB,GPX3,  SRC 

hypertrophy  of  cardiomyocytes 

6.1 6E  06 

A^1,EP300,HDAC2,LEP,MAPK1,PI3K  (complex),  STAT3 

coronary  artery  disease 

3.08E  04 

APOE,  CNR1,E0S,HDAC2,KL,NR3C1,  LEP 

metabolic  disorders:  types  1  and  II  diabetes 

3.78E  15 

AANAT,ABCC8,ACLY,ACSBG1  ,ACSBG2,ACSL4,ACSL6,ACVR1  C,ADCY8,ADIPORl , 
ADIPOR2,AGRP,A^l  ,ANXA1  ,APAE1  ,APOE,ASIC2,ATM,BAD,BCARl  ,CACNA1  C, 
CAMP,CAPN8,CASP3,CBL,CCKAR,CD247,CDK2,CEBPB,CERK,CHUK,CITED1,CNR1, 
CNR2,CPE,CPLX3,GAD2,GCK,GH,GHRHR,GHRL,GHSR,GLP1  R,GPR1 1 9,GRB2,GSK3B, 
GYSl  ,GZMB,H2DMA,H2DMB1  ,H2Q1  ,H2BEM,HADH,HNE1  A,HPSE,HSPD1  ,H^,ICA1 , 
IENGR2,IKBKG,ILDR2,INPP5K,INPPL1,INSR,IRAK1,IRS3,JAK2,KCNJ11,KHK,KRAS,LARS, 
LEP,LEPR,LEPR0T,LIPE,LMAN1L,LPIN1,LTA,MAP2K2,MAP3K7,MAPK8,MC4R,MRAS, 
MT0R,MY05A,NCK1  ,NGER,NOS2,NPY,NRAS,NSMAE,OCRL,PAKl  ,SGK1  ,SH2B2,SIK2, 
SLC2A4,SLC2A8,SLC30A8,SMAD2,SOCS1  ,SORBSl  ,SORTl  ,SOSl  ,SOS2,SRC,STATl , 
STAT3,STX1  A,STXBP4,SYNJ2,SYP,TLN  1  ,TLN2,TRADD,TRAE6,TRIB3,TRIP1 0,TSC1 , 
TSC2,VAMP2,VDR,WDTC1  ,ZEP1 06 

weight  gain  and  obesity 

4.33E  10 

APP,BDNE,CNR1,CRH,DRD2,E0S,LEP,NPY,NTRK2,  STAT3,CRK,CYB5R4,CYP27B1, 
DGATl  ,DOK4,EGRl  ,EIE2AK3,EIE4EBP2,ENPP7,EADD,EAM3B,EAS,EASL,ECER1  G, 
EEAR1,EKBP1B,E0XC2,E0X04,ERS3,EYN,GAB1,PCL0,PCSK9,PDE3B,PDIA3,PDK1, 
PDPKl  ,PDX1  ,PHIP,PIAS1  ,PIK3C2G,PIK3CG,PIK3R5,PKLR,PKM,PLCB2,PLCZ1  ,PPARG, 
PPP1R3D,PRE1,PRKAR2B,PRKCD,PTEN,PTPN2,PTPRA,RAE1,RAPGEE4,RBP4,RELA, 
RETN,RH0Q,RIMS2,RIPK1,RPE65,RPS6KB2,RPT0R,RRAS2,SCAP,SCNN1G, 

set  of  DEGs,  the  expression  patterns  of  AY  and  HC  were 
reversed  for  the  early  (5dld)  and  for  the  late  (10d42d) 
time  points  for  transcripts  filtered  with  stricter  cut  offs. 


This  corroborates  our  suggested  activation  and  inhib¬ 
ition  of  behavioral  responses  across  time  points  among 
AY,  HC  and  MPFC.  Transcripts  at  early  time  point 
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Figure  8  Transcripts  brain  regions  interactome  for  transcripts  that  passed  p  value  less  than  0.01  and  fold  change  greater  than  2.0  in 
Agg  E  groups  compared  to  control  groups  (across  time  points). 


(5dld)  from  AY  were  largely  induced  whereas  they  were 
suppressed  in  HC  (Figure  8,  Additional  file  1:  Figure  S9). 
On  the  other  hand,  transcripts  at  later  time  point 
(10d42h)  showed  reversed  expression  pattern  between  AY 
and  HC  (Figure  8,  Additional  file  1:  Figure  S12),  which  is 
again  consistent  with  our  expression  based  prediction  ana¬ 
lyses.  At  the  transcript  level,  there  were  not  many  shared 
DEGs  among  combinations  of  brain-regions  and  time 
points.  But  there  were  more  overlaps  at  the  pathway  and 
process  levels  (Additional  file  1:  Figures  S9-S18). 

A  large  portions  of  these  sets  of  transcripts  were  asso¬ 
ciated  with  splice  variant  processing  and  activities  of 
zinc  finger  transcription  factors  (Figure  9  and  Additional 
file  3:  Table  S2).  Many  of  these  transcripts  were  also 
shared  also  with  those  transcripts  suggested  to  be  in¬ 
volved  in  traumatic  disorders  and  comorbid  conditions 
(Figure  9  and  Additional  file  3:  Table  S2). 

Validation  of  gene  expression  profiles 

Transcripts  implicated  in  anxiety  and  fear  responses  in 
PTSD  patients  and  animal  models  were  assayed  in  each 
brain  region  (AY,  HC,  MPFC,  SE,  ST  and  VS)  at  the 
TlORl  time  point  using  the  QuantiGene  Plex  2.0  Reagent 
System  (Multiplex)  to  corroborate  our  microarray  results. 


Many  of  the  samples  assayed  by  the  QuantiGene  Multi¬ 
plex  showed  comparative  direction  of  expression  with 
that  of  microarray  data  (Additional  file  1:  Figures  S19.a 
and  S20).  Transcripts  critical  in  nerve  growth  factor  and 
synaptic  plasticity  such  as  FK506  binding  protein  5 
(Fkbp5),  brain-derived  neurotrophic  factor  (Bdnf),  thy¬ 
roid  hormone  receptor  beta  (Thrb)  and  prostaglandin- 
endoperoxide  synthase  2  (Ptgs2)  were  largely  suppressed 
in  both  microarray  and  QuantiGene  assays.  Transcripts 
of  the  neuronal  signaling  molecules  (substrate  specific 
channels)  such  s  Slcla2  was  up-regulated  and  Ill  rap  was 
down-regulated  in  all  regions  except  in  AY,  while  Rx^3 
and  Proxl  were  up-regulated  in  AY,  HC  and  MPFC  re¬ 
gions  (Additional  file  1:  Figures  S19.a  and  S20). 

Correlation  analysis  for  each  transcript  between  the 
bead  and  array  data  showed  directional  correlation  for 
the  transcripts:  Fkbp5,  Thrb,  Illrap,  Bdnf,  Drd2,  Slcla2, 
Coch,  Gng4,  Rxfp3  and  Slcla2;  while  Gpx3,  Npy,  Proxl 
and  Ptgs2  showed  lower  directional  correlation  between 
the  two  platforms  (Additional  file  1:  Figures  S19  and  S20). 
Rxfp3,  Gng4  and  Slcla2  showed  overall  opposite  correl¬ 
ation  compared  to  the  rest  of  the  transcripts  as  shown 
by  the  correlation  matrix  of  all  the  assayed  transcripts 
(Additional  file  1:  Figure  S19.b).  These  transcripts  were 
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Figure  9  Pathways  and  biological  processes  significantly  associated  with  transcripts  (with  p  value  <  0.01,  and  fold  change  >  2.0)  shown 
in  Figure  8.  Sizes  of  the  hollow  pie  sectors  were  proportional  to  the  number  of  transcripts  associated  with  that  particular  pathway  or  process. 
Gene  expression  related  processes  such  as  activities  of  zinc  finger  protein  transcription  factors,  and  splice  variant  processing  were  associated  with 
larger  numbers  of  transcripts  (and  many  of  these  were  also  shared  with  pathways  related  to  trauma/fear  related  behavioral  responses). 


mostly  up-regulated  in  both  platforms  (red  ellipses) 
where  as  the  other  transcripts  showed  similar  correl¬ 
ation  and  were  largely  down-regulated  in  the  two  plat¬ 
forms  (blue  ellipses)  (Additional  file  1:  Figure  S19.b). 
This  also  corroborated  the  notion  that  neurogenesis  re¬ 
lated  brain  factors  were  suppressed  whereas  some  of 
transcripts  implicated  substrate  specific  channels  and 
neuronal  signaling  were  induced  at  the  molecular  level. 

Discussion 

Cellular  and  molecular  changes  in  the  activity  and  func¬ 
tional  connectivities  among  the  amygdala  (AY),  hippo¬ 
campus  (HC),  medial  prefrontal  cortex  (MPFC)  play 
primary  roles  in  fear  learning  and  memory,  and  thus. 


may  contribute  to  the  behavioral  manifestations  of  anx¬ 
iety  disorders  (Additional  file  1:  Figure  SI). 

Identification  of  the  molecular  basis  of  learning  and 
memory  mechanisms  across  stress-induced  anxiety  disor¬ 
ders  is  critical  for  devising  means  of  mitigating  the  severity 
of  traumatic  disorders.  This  is  because  both  memory  re¬ 
trieval  (involved  in  symptom  expression,  reconsolidation, 
and  maintenance  of  these  disorders),  and  memory  extinc¬ 
tion  (believed  to  be  the  mechanism  of  behavioral  exposure 
therapy  of  anxiety  disorders)  are  dependent  on  transcrip¬ 
tional  changes  of  the  underlying  molecules  governing  syn¬ 
aptic  plasticity  and  inter-neuronal  signaling. 

Transcripts  significantly  associated  with  synaptic  transmis¬ 
sion  and  plasticity,  long-term  depression  and  potentiation. 
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axonogenesis,  and  the  synthesis  and  release  of  neurotrans- 
mitters  such  as  glutamate,  acetylcholine,  dopamine  and 
serotonin  were  differentially  regulated  across  time  points 
and  brain  regions  due  to  the  effects  of  social  stress  trauma 
on  Agg-E  mice.  These  molecular  changes  are  corroborated 
by  changes  in  other  groups  of  transcripts  associated  with 
behavioral  expression  of  anxiety  disorders  such  as  startle 
response,  anxiety,  hyperactivity,  fear  conditioning,  emo¬ 
tional  behavior,  cognitive  impairment,  and  addiction  and 
alcohol  abuse-associated  behaviors.  Molecular  indicators 
of  addiction  and  alcohol  abuse  showed  increased  activa¬ 
tion  z-scores  (Additional  file  1:  Figure  S7A),  which  cor¬ 
relates  with  the  observations  of  others  that  traumatic 
events  sufficient  to  produce  long-lasting  enhancement 
of  fear  learning  also  increase  voluntary  ethanol  con¬ 
sumption  [36]. 

DEGs  involved  in  synaptic  transmission  in  relation  to 
long-term  synaptic  potentiation,  long-term  synaptic  de¬ 
pression,  and  synaptic  plasticity  are  the  basis  for  learning 
and  memory  consolidation.  These  pathways  and  func¬ 
tional  enrichment  findings  may  indicate  the  acquisition, 
consolidation  and  maintenance  of  traumatic  memories 
to  form  anxiety  disorders  via  fear  learning  and  fear  mem¬ 
ory  circuitry,  as  well  as  deficiencies  in  extinction  and 
re-extinction  mechanisms  along  the  different  time  course 
trajectories  in  different  brain  regions.  For  example,  the  AY 
had  higher  activation  z-scores  of  anxiety- related  behaviors 
at  earlier  time  points  compared  to  the  MPEG,  which  cor¬ 
responded  with  the  earlier  response  in  the  amygdala.  This 
supports  the  belief  that  neural  circuits  underlying  PTSD 
are  characterized  by  exaggerated  AY  activation  and  de¬ 
creased  MPEG  activation,  leading  to  the  elevated  anxiety 
state  and  concomitant  inadequate  emotional  regulation 
[37].  Re-experiencing  the  traumatic  event  in  the  form  of 
recurrent  distressing  images  and  recollections,  including 
the  intrusive  traumatic  memory,  is  unique  to  PTSD 
among  anxiety  disorders  [17,38],  and  is  probably  due  to 
diversion  of  focus  from  the  current  stimuli  (caused  by  im¬ 
paired  cognition  due  to  slow  acting  MPEG)  to  the  intru¬ 
sive  traumatic  event  (interference  effect  of  the  faster  paced 
AY  activities). 

In  addition,  DEGs  reported  to  be  involved  in  the  fear- 
potentiated  startle  response,  fear  extinction  resistance, 
associative  fear  learning,  and  cue  conditioning  were  indi¬ 
cative  of  the  formation  of  traumatic  memories  and  a  vul¬ 
nerability  to  developing  fear-related  anxiety  disorders.  For 
example,  the  transcripts  FKBP5  and  BDNF  were  sup¬ 
pressed  across  the  six  brain  regions  at  TlORl  (Additional 
file  1:  Figures  S19.a  and  S20).  FKBP5  is  involved  in  the 
behavioral  and  neuroendocrine  effects  of  chronic  social 
defeat  stress.  Suppression  of  FKBP5  expression  in  HG 
leads  to  persistent  traumatic  disorder  and  confers  fear  ex¬ 
tinction  resistance  [22,23].  BDNF  is  a  crucial  regulator  of 
neuroplasticity,  which  underlies  learning  and  memory 


processes  in  different  brain  areas.  Ghronic  suppression  of 
BDNF  expression  in  the  HG,  AY  and  the  prefrontal  cortex 
leads  to  a  deficit  in  the  acquisition  of  extinction  memory, 
and  increased  BDNF  expression  supported  amelioration  of 
hyperarousal  in  a  mouse  model  of  PTSD  [39].  BDNF  facil¬ 
itates  extinction  learning,  and  decreased  mRNA  expres¬ 
sion  of  BDNF  within  the  MPFG  leads  to  resistance  to  fear 
extinction.  Activation  of  BDNF  signaling  by  the  TrkB 
agonist  7,  8-dihydroxyflavone  blocks  the  return  of  fear 
after  extinction  training  [40,41]. 

The  exact  social  stress  process  that  we  used  here  was 
also  shown  to  result  in  decreased  BDNF  in  the  HG,  which 
led  to  suboptimal  binding  of  BDNF  to  tyrosine  receptor 
kinase  B  (TrkB),  resulting  in  curtailing  downstream  intra¬ 
cellular  signaling  pathways,  including  mitogen-activated 
protein  kinase/extracellular  signal- regulated  protein  kinase 
(MAPK/ERK),  phospholipase  Gg  (PLGg),  and  phosphoino- 
sitide  3-kinase  (PI3K)  pathways.  These  pathways  are  im¬ 
portant  for  neurogenesis  and  cognition  [42,43].  Another 
important  DEG  along  this  pathway  is  the  nerve  growth 
factor,  Ngf,  that  also  shown  to  bind  to  tyrosine  receptor 
kinase  A  (TrkA).  In  response  to  Ngf,  the  TrkA  forms  a 
complex  with  She,  coupling  TrkA  to  p21ras  and  She  with 
Grb2,  which  is  mediated  by  autophosphorylation  of  TrkA 
[44].  However,  suppression  of  Ngf  at  earlier  time  points 
(T5R1  and  TlORl)  suggestive  of  curtailed  kinase  activities 
and  complex  formations,  potentially  leading  to  impaired 
synaptogenesis. 

Potential  contribution  to  the  etiology  of  PTSD  may  be 
due  to  sensitization  of  glucocorticoid  receptor  (GR)  sig¬ 
naling  and  dysregulation  of  GR  modulators  [45]  such  as 
FKBP5  and  Grfl  receptors.  Furthermore,  Akt,  Nfkb  and 
MAP  kinases,  which  are  G  protein-coupled  receptor 
(GPGR)  pathway  molecules,  can  promote  or  prevent 
sustained  high  anxiety-  and  depressive-like  behaviors  fol¬ 
lowing  severe  stress.  Agonist-induced  activation  of  the 
corticotropin  releasing  factor  (Grfl  receptor)  is  crucial 
for  survival  in  the  context  of  serious  danger  or  trauma, 
but  persistent  Grfl  receptor  hyper- signaling  when  a 
threatening  or  traumatic  situation  is  no  longer  present  is 
maladaptive.  Also,  suppressed  expression  of  Grk3,  which 
phosphorylates  the  Grfl  receptor  protein  indicates  there 
was  suboptimal  binding  of  Grhrl  with  beta-arrestin2, 
leading  to  inhibited  termination  of  Grk3-coupled  Grfl 
receptor  signaling  by  homologous  desensitization.  This 
deficiency  in  Grk3-beta-arrestin2  complex  formation  is 
suggested  to  contribute  to  PTSD  and  co-morbid  post- 
traumatic  depression  [46]. 

FKBP5  is  also  shown  to  mitigate  HPA  axis  dysfunction 
[47].  Molecular  changes  involved  in  modulation  of  the 
HPA  axis  are  associated  with  cue  conditioning,  calcium 
signaling,  memory  processes,  and  regulation  of  the  stress 
responses.  Transcripts  of  genes  involved  in  HPA  axis, 
such  as  those  related  to  corticotrophin  releasing  factor 
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(CRF)  receptors  pathways,  cortisol  (glucocorticoid)  re¬ 
ceptor  signaling  and  oxytocin  receptor  (such  as  NrScl, 
Nr3c2,  Crhrl,  Crhr2,  Oxtr)  were  also  differentially  regu¬ 
lated.  Oxytocin  (mammalian  neuropeptide)  modulates 
activation  of  fear  extinction-based  neural  circuits  and 
fear  responses  [48]. 

Release  of  GC,  CRF  and  modulators  are  shown  to 
affect  emotional  learning  and  memory  in  relation  to  eti¬ 
ologies  of  anxiety  disorders.  The  pathogenesis  of  PTSD 
is  attributed  to  over-consolidated  traumatic  memories 
that  are  mediated  by  endogenous  stress  hormones  re¬ 
leased  during  and  after  trauma.  Downstream  secondary 
messenger  signaling  pathways  with  a  putative  role  in 
long-term  potentiation  (LTP),  such  as  inositol  1,4,5-tris- 
phosphate  (IPS)  and  diacylglycerol  (DAG)  with  the  up¬ 
stream  enzyme  inositol  monophosphatase  (Imp),  favor 
fear  learning  and  traumatic  memory  consolidation  in 
traumatized  and  susceptible  individuals  via  fear  potenti¬ 
ated  LTP,  and  the  eventual  development  of  PTSD  [49]. 
Glucocorticoid  increase  after  longer  rest  periods  in  Agg- 
E  mice  is  probably  an  indicator  of  attenuated  condi¬ 
tioned  fear  and  is  in  agreement  with  the  observation  that 
stress-induced  circulating  cortisol  was  shown  to  reduce 
memory  retrieval  of  conditioned  fear  in  men  [50]. 

Transcript  levels  of  genes  reported  to  be  associated 
with  cognition  and  object  recognition  (in  particular,  sup¬ 
pressed  expression  of  BDNF,  Gnail,  App,  EpSOO,  Cnrl, 
Cnr2  and  Ngf  transcripts)  at  one-day  post- 10  days  of 
Agg-E  session  (TlORl)  indicated  there  was  impaired 
cognition,  and  therefore,  inhibition  of  novel  object  rec¬ 
ognition.  The  memory  modulating  cannabinoid  recep¬ 
tors  (Cnrl  &:  Cnr2)  and  elevated  levels  of  endocannabinoids 
(endogenous  cannabinoids)  in  the  basolateral  AY  and  have  a 
critical  function  in  the  extinction  of  aversive  memories. 
Cnrl -deficient  mice  showed  strongly  impaired  short-term 
and  long-term  extinction  in  auditory  fear-conditioning 
tests,  with  unaffected  memory  acquisition  and  consolida¬ 
tion.  Treatment  of  wild-type  mice  with  the  Cnrl  (CBl 
receptor)  antagonist  SR141716A  mimicked  the  phenotype 
of  CBl -deficient  mice,  suggesting  that  CBl  is  required 
at  the  moment  of  memory  extinction.  In  the  basolateral 
AY,  endocannabinoids  and  CBl  were  crucially  involved 
in  long-term  depression  of  GABA  (gamma-aminobutyric 
acid) -mediated  inhibitory  currents  [51,52]. 

Alternatively,  transcripts  of  the  histone  deacetylases 
(Hdac2  &:  HdacS)  showed  increased  expression,  also  indi¬ 
cating  impaired  object  recognition.  Impaired  object  recog¬ 
nition  memory  may  underlie  certain  avoidant  symptoms 
or  negative  cognitions  in  PTSD  and  be  related  to  impaired 
behavioral  flexibility  [53].  Impaired  novel  object  recogni¬ 
tion  parallels  neurocognitive  deficit  (impaired  memory 
and  attention)  in  PTSD  patients  who  show  lower  perform¬ 
ance  on  memory  tests  specific  to  learning  and  executive 
functions  [54]. 


Other  important  DEGs  include  cFos  and  its  (pro¬ 
moter)  histone  deacetylase  enzyme,  Hdac2,  GluN2B 
(Duspl),  and  klotho  (/<7).  Hdac2  and  cFos  play  a  critical 
role  in  fear  memory  recall  for  reconsolidation  and  up¬ 
dating  acquired  memories  in  neuronal  plasticity  [38,39]. 
A  similar  social  defeat  model  showed  that  presentation 
of  an  aggressor  cue  induced  robust  increases  in  cFos  in 
AY,  MPFC  and  CAl  of  HC  both  1  and  7  days  after  Agg- 
E.  The  increase  at  7  days  is  greater  than  that  at  one  day, 
whereas  avoidance  behavior  was  great  at  one  day  post- 
Agg-E  and  had  somewhat  abated  by  7  days  post- Agg-E. 
This  suggests  that  increased  neuronal  processing  was  re¬ 
lated  to  the  decrease  in  avoidance  [55].  In  the  condi¬ 
tioned  fear  paradigm  of  mice,  fear  memory  extinction  is 
shown  to  be  more  effective  in  a  one-day-old  fear  mem¬ 
ory  compared  to  a  30-day-oId  fear  memory.  Older  fear 
memories  are  less  labile  due  to  hypoacetylated  c-Fos 
promoter  by  Hdac2  in  HC  [56].  Among  others,  this  may 
be  one  reason  why  we  observed  more  DEGs  after  42  days 
of  rest  (T10R42)  compared  to  one-day  post  10-day  Agg- 
E  sessions  (TlORl),  with  related  DEGs  activated  upon 
remote  memory  recall  [39,40].  Elevated  K1  in  mice  has 
been  shown  to  enhance  long-term  potentiation  (synaptic 
plasticity),  and  enrich  synaptic  GluN2B  (an  NMDA  re¬ 
ceptor  subunit)  with  key  functions  in  learning  and  mem¬ 
ory.  Blockade  of  GluN2B  abolished  klotho-mediated 
effects.  Suppressed  expression  of  K1  may  indicate  cogni¬ 
tive  deficits  in  our  model  animals  [57]. 

Overall,  transcripts  of  DEGs  in  Agg-E  groups  of  mice 
compared  to  controls  in  the  different  brain  regions  were 
reported  to  be  associated  with  glucocorticoid  negative 
feedback  signaling,  arousal  to  trauma  cues  (jumping),  im¬ 
paired  aggressive  behavior  (avoidance  of  aggressor-cued 
partition),  social  withdrawal,  impaired  territorial  behavior, 
long-term  fear  memory,  grooming,  decreased  movement, 
anxiety,  long-term  synaptic  depression  and  potentiation 
(long-term  synaptic  plasticity),  inhibited  dopaminergic 
signaling,  HC  atrophy,  dendritic  branching  in  AY,  condi¬ 
tioned  fear  association,  fear  memory  consolidation,  retrieval 
and  impaired  contextualization,  sensorimotor  gating  deficit, 
deficit  of  executive  function,  impaired  object  recognition, 
and  circadian  rhythm  disruption  (Table  2). 

Significant  association  of  DEGs  with  comorbid  condi¬ 
tions,  such  as  chronic  inflammation,  myocardial  infarc¬ 
tion,  suppressed  protective  immunity  and  obesity/ 
diabetes  indicate  the  pervasive  nature  of  Agg-E,  leading 
to  many  systemic  pathological  consequences.  Previous 
studies  using  this  same  mouse  model  simulating  PTSD 
showed  acute  myocardial  injury  associated  with  the  trau¬ 
matic  experience  as  a  consequence  of  underlying  bio¬ 
logical  injury  processes,  including  inflammation  [58]. 
Metabolomic,  histopathology  and  liver  transcriptomic 
analyses  showed  increased  inflammatory  response  at 
one-day  post-stress,  persistent  myocarditis  and  cardiac 
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fibrosis,  and  hyperlipidemia  in  the  liver,  related  to  risks 
of  heart  disease  and  obesity  [59].  Similar  co-morbidities 
of  PTSD  observed  in  veterans  and  active  duty  personnel 
include  inflammation,  obesity,  diabetes,  and  heart  dis¬ 
ease  [3],  as  well  as  sleep  disorders  [60,61]  and  higher 
rates  of  pain,  possibly  due  to  chronic  inflammation  [62] . 

Altered  trauma  induced  behavioral  responses  and  as¬ 
sociated  comorbid  pathways  (discussed  before  and  add¬ 
itional  pathways)  were  also  significantly  associated  with 
the  more  significant  transcripts  (with  p  <  0.01,  and  fold 
change  >  2.0)  (Figure  9  and  Additional  file  1:  Figures  S9- 
S18).  Additional  pathways  related  to  splice  variant  pro¬ 
cessing,  mitochondrial  dysfunction,  metabolic  disorder, 
substrate  specific  channel  activities,  activities  of  zinc  fin¬ 
ger  related  transcription  factors  and  cytoskeletal  protein 
binding  were  significantly  enriched  (Figure  9).  Associ¬ 
ation  of  many  of  the  top  transcripts  with  splice  variant 
processing  may  hint  to  the  basis  of  molecular  mecha¬ 
nisms  leading  to  either  stress  resilience  or  susceptibility 
among  individuals.  Transcripts  associated  with  mito¬ 
chondrial  dysfunction  along  with  those  involved  in 
metabolic  disorder  may  contribute  towards  exaggerated 
fear  in  traumatized  individuals  [63]. 

Conclusions 

Pathways  and  functions  presumably  important  for  the 
etiology  of  PTSD-like  disorders,  including  endocannabi- 
noid  signaling,  HPA  axis  function,  modulators  and  tar¬ 
gets  of  cortisol,  neuronal  transmission,  neurogenesis  and 
fear  memory  extinction  with  regard  to  emotional  learn¬ 
ing  and  memory,  and  a  number  of  anxiety-related  behav¬ 
ioral  responses  were  significantly  associated  with  DEGs. 

Activation  z-scores  of  anxiety- related  behavioral  re¬ 
sponses,  fear  acquisition  and  consolidation-related  sig¬ 
naling  pathways  and  processes  reveal  the  distinct,  and 
probably  complementary  roles  of  early  time  points  ver¬ 
sus  later  rest  periods  (of  5-day  and  10-day  Agg-E  ses¬ 
sions)  in  the  traumatic  fear-learning  trajectories.  The 
molecular  mediators  of  the  earlier  stress  responses  (first- 
wave  responders)  seem  to  enhance  fear  maintenance  (in 
AY),  whereas  the  second-wave  of  transcripts  seem  to  ei¬ 
ther  consolidate  fear  memory  (in  HC)  or  attenuate  fear 
learning  (in  MPFC). 

Future  directions  and  suggested  works 

Advanced  experimental  tools,  such  as  inducible  muta¬ 
tions  in  mice,  virus -mediated  gene  transfer,  and  optoge- 
netics  provide  the  means  to  directly  delineate  the  role  of 
specific  proteins  acting  within  specific  cell  types  of 
anxiety-processing  brain  structures  in  mediating  PTSD- 
like  behavioral  abnormalities  in  animal  models. 

There  is  evidence  that  females  use  different  cognitive 
strategies,  exhibit  increased  stress  sensitivity,  and  show 
variations  in  reward-related  behaviors  throughout  the 


estrus  cycle  that  may  render  them  more  sensitive  to  the 
deleterious  effects  of  stress.  Future  studies  may  examine 
the  molecular  and  cellular  underpinnings  of  anxiety-like 
behaviors  in  females.  For  example,  depression  is  twice 
as  likely  to  occur  in  women  as  in  men  [24],  and  this 
may  be  true  for  other  anxiety  disorders  as  well,  but  ani¬ 
mal  studies  including  ours  have  mostly  been  conducted 
in  males. 

Untreated  traumatic  fear  memory  in  susceptible  indi¬ 
viduals  precipitates  in  the  form  of  anxiety  disorders  such 
as  PTSD.  One  effective  treatment  of  anxiety  disorder  is 
using  exposure  therapies,  which  are  essentially  traumatic 
fear  memory  extinction  and  are  effective  during  the  la¬ 
bile  phase  of  fear  memory  consolidation.  Since  the  win¬ 
dow  of  fear  memory  consolidation  is  relatively  short 
compared  to  the  enduring  flashback  vulnerability,  re¬ 
mote  fear  memories  may  not  be  persistently  unlearned 
by  reconsolidation-updating  paradigms.  Hence,  there  is 
a  need  to  identify  molecular  mediators  of  stable  fear 
memory  in  order  to  devise  a  mechanism  of  extinction  of 
remote  fear  memory  long  after  its  consolidation  window. 

Methods 

Exposure  of  subject  mice  to  aggressor  mice 

Albino  SJL  male  mice  (The  Jackson  Laboratory,  Bar 
Harbor,  ME)  were  single-housed  for  one  month  and 
trained  as  aggressors  by  placing  bulbectomized  C57BL/6 
male  mice  in  SJL  cages  as  intruders.  [Note:  the  purpose 
of  bulbectomized  C57BL/6J  mice  (which  don’t  fight 
back)  was  only  to  train  the  SJL  mice  to  be  more  ag¬ 
gressive].  Naive  C57BL/6J  male  mice  single-housed  in  a 
separate  room  were  used  as  experimental  subjects  (5 
exposed  and  5  control  mice  per  time  point  were  used). 
Each  subject  mouse  was  kept  in  a  protective  mesh  box 
inside  the  cage  of  the  aggressor  SJL  mouse  for  6  hours/ 
day  for  5  (T5)  or  10  days  (TIO),  with  direct  exposure  to 
the  aggressor  mouse  for  either  1  minute  or  10  bites 
(whichever  occurs  first)  for  1-3  random  times  during 
the  6-hour  defeat  session.  Control  (C-ctrl)  mice  were 
also  housed  without  feed  and  water  for  the  same  period 
in  a  room  separate  from  the  aggressor  mice. 

Body  weights  and  temperatures  of  the  experimental 
mice  were  measured  (using  implanted  electronic  ID 
transponders)  before  and  after  the  social  defeat  session 
of  each  day.  Details  are  given  in  separate  publications 
[29].  Territorial  urine  markings  were  measured  by  pla¬ 
cing  blotting  papers  (0.8  mm  thick)  under  the  cages  of 
Agg-E  and  C-ctrl  mice.  Areas  were  compared  as  de¬ 
scribed  in  [29]. 

Assays  for  corticosterone  levels 

One  hour  after  the  partition  test,  mice  were  euthanized  by 
cervical  dislocation  and  blood  samples  were  immediately 
drawn  by  cardiac  puncture  and  collected  in  nuclease-free 
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tubes  containing  50  \i\  of  buffered  sodium  citrate 
(0.105  M  ~  3.2%).  Whole  blood  samples  were  then  centri¬ 
fuged  and  plasma  was  separated  and  stored  at  -80°C  until 
the  corticosterone  assay  was  performed.  Corticosterone 
level  was  measured  using  an  ELISA  kit  following  the  ma- 
nufacturers  instructions  (Kamiya  Biomedical  Company, 
Seattle,  WA). 

Brain  dissection  and  collection  of  brain  regions 

Brain  regions  from  subject  and  control  C57BL/6  mice 
were  collected  after  one  (Rl)  or  42  days  (R42)  of  post 
10-day  (TIO)  aggressor  exposure,  and  one  (Rl)  or  10  days 
(RIO)  after  the  5-day  (T5)  aggressor  exposure.  Brains 
were  carefully  removed  from  the  skulls,  and  the  left  or 
right  hemi-brain  from  each  traumatized  or  control 
mouse  was  dissected  into  different  anatomical  and  func¬ 
tional  regions:  amygdala  (AY),  hippocampus  (HC),  med¬ 
ial  prefrontal  cortex  (MPFC),  ventral  striatum  [including 
nucleus  accumbens]  (VS),  septal  region  (SE)  and  corpus 
striatum  (ST).  Collected  tissues  were  immediately  frozen 
in  liquid  nitrogen  and  stored  at  -80°C  for  total  RNA 
isolation. 

RNA  isolation 

Total  RNAs  were  isolated  according  to  the  Trizol®  method 
(Invitrogen,  Inc.,  Grand  Island,  NY)  from  homogenized 
brain  regions.  Isolated  RNAs  were  stored  at  -80°C  until 
microarray  assays  were  conducted. 

Microarray  hybridization 

Microarray  assays  were  performed  using  Agilents 
genome-wide  mouse  expression  array  (GE  4x44K  v2  two 
color  microarray)  slides  and  kits  (Agilent  Technologies, 
Inc.,  Santa  Clara,  CA)  following  the  manufacturer  s  proto¬ 
col.  Hybridized  microarray  slides  were  scanned  using 
Agilent  Technologies  Scanner  G2505C  US09493743. 

Feature  extraction  and  normalization 

Images  of  scanned  microarray  slides  were  feature- 
extracted  and  normalized  using  Agilent  s  feature  extrac¬ 
tion  software,  version  10.7  or  later,  in  the  default  setup 
(Agilent  Technologies,  Inc.). 

Microarray  data  analysis 

Normalized  data  were  filtered  on  flags  to  exclude  probes, 
which  have  missing  values  in  more  than  one  sample,  and 
these  were  tissue-wise  quantile  normalized  using  the 
Limma  Package  of  R  program  (www.r-project.org).  Nor¬ 
malized  data  were  analyzed  using  Agilents  Genespring 
GX  vl2.0,  the  Limma  package  of  R,  and  Qlucore  Omics 
Explorer  2.3  (Qlucore  AB,  Lund,  Sweden).  Each  brain  re¬ 
gion  at  each  time  point  (T5R1,  T5R10,  TlORl,  and 
T10R42)  was  compared  for  aggressor  exposure  effects: 
Aggressor-Exposed  (Agg-E)  vs.  Control  (C-ctrl).  DEGs 


in  each  tissue  and  time  point  were  identified  using  a 
Moderated  T-test  2it  p  <  0.05  of  the  Limma  Package  of  R. 

Time  effects  were  also  compared  using  TimeClust 
[64]  across  different  combinations  of  aggressor  expos¬ 
ure  sessions  and  rest  periods  (T5R1,  T5R10,  TlORl,  and 
T10R42). 

Behavioral  association  of  DEGs  with  respect  to  functional 
and  pathway  enrichments 

To  assess  behavioral  association  of  DEGs,  we  followed 
two  complementary  approaches: 

i.  We  started  with  DEGs  across  different  brain  regions 
and  identified  significantly  enriched  signaling 
pathways  and  biological  processes. 

ii.  We  used  frequency  calculation  to  determine 
psychopatholy-related  pathways  and  biological 
processes  cited  in  the  scientific  literature  in  relation 
to  social  stress,  trauma,  PTSD,  anxiety,  and  depression. 
We  then  identified  (from  databases  and  literature) 
transcripts  and  proteins  directly  associated  with 
the  most  frequently  cited  terms.  Literature-  and 
database-mined  genes  and  proteins  were  intersected 
with  DEGs  from  our  study  (brain  regions),  and 
those  in  common  were  used  to  determine  probable 
downstream  biological  processes,  as  well  as  networks 
and  pathways  and  their  significant  associations. 

Gene  ontology  and  pathway  enrichments,  regulatory 
networks  analysis,  common  pathways  analysis,  and 
visualization 

The  Hypergeometric  Test  (false  discovery  rate  (FDR), 
q  <  0.05)  of  Bingo  2.44  and  ClueGO  (Cytoscape  3.0.1  plug¬ 
ins;  http://www.cytoscape.org/),  Fishers  Exact  Test  of 
Ingenuity®  (Ingenuity,  Inc.,  Redwood,  CA),  and  Gene  Set 
Enrichment  Analysis  (http://www.broadinstitute.org/gsea/ 
index.jsp)  were  used  for  gene  ontology  and  pathway  en¬ 
richments.  The  Cytoscape  and  Ingenuity  Pathway  analyt¬ 
ical  software  systems  were  also  used  for  network  analyses 
and  visualizations.  Common  pathways  were  intersected 
and  visualized  using  Cytoscape  and  graphically  displayed 
using  the  R  programming  (www.cran-r.org)  graphic  system. 
Activation  z-scores  of  behavioral  responses  and  synaptic 
plasticity  were  used  to  identify  induced  and  suppressed 
behaviors  and  synaptic  activities  across  brain  regions 
and  time  points. 

Transcript  -  brain-region  interactomes  were  constructed 
and  visualized  using  Gephi.0.8.2  beta  (www.gephi.org). 
Functions  and  pathways  associated  with  the  top  transcripts 
(p-value  <  0.01,  and  fold  change  >  2.0)  were  enriched  using 
ClueGo,  Bingo,  IPA  and  David  (using  FDR  correction  of 
q  <  0.05).  Bar  charts  of  the  corresponding  processes  and 
pathways  were  constructed  using  the  ggplot2  of  R. 
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Gene  expression  profiling  using  QuantiGene  Plex  2.0 
Reagent  system 

Target- specific  RNA  molecules  {Thrb,  NM  009380;  Ptgs2, 
NM  011198;  Proxh  NM  008937;  Drd2,  NM  010077; 
Slcla2,  NM  011393;  Rgs4,  NM  009062;  Hmga2,  NM 
010441;  Npy,  NM  023456;  Rxp3,  NM  178717;  Gng4, 
NM  010317;  NM  007540;  Coch,  NM  007728;  Illrap, 
NM  134103;  NM  008161;  FkbpS,  NM  010220) 

from  different  brain  regions  were  analyzed  by  the  Quanti¬ 
Gene  Plex  2.0  Reagent  system  (Affymetrix,  Fremont,  CA). 
Oligonucleotide  probe  sets  (capture,  label,  and  blocker 
probes)  for  each  gene  were  designed  by  the  manufacturer. 
Briefly,  RNAs  from  different  brain  regions  were  captured 
by  fluorescent  microspheres.  After  overnight  hybridization 
at  54  ±  1°C,  hybridizations  with  branched  DNA  (bDNA) 
pre-amplifier  2.0,  bDNA  amplifier  2.0,  biotinylated  label 
probe  and,  finally  substrate  were  sequentially  carried  out 
according  to  the  manufacturer  s  instructions.  Signals  of  cas¬ 
cade  amplification  were  detected  by  Bio-Plex  100  x  MAP 
technology  and  analyzed  using  Bio-Plex  6.0  software 
(Bio-Rad  Laboratories,  Hercules,  CA).  Mean  fluorescence 
intensity  (MFI)  signals  generated  from  each  bead  are  pro¬ 
portional  to  the  amount  of  each  mRNA  captured  on  the 
surface  of  each  generated  specific  probe  set  [65]. 

The  geometric  means  of  four  housekeeping  genes  (Ppib, 
NM  011149;  HPRTl,  NM  013556;  Ldha,  NM  010699; 
RplpO,  NM007475)  were  used  for  normalization  of  each 
sample.  Fold-changes  were  the  relative  ratios  between  nor¬ 
malized  values  of  treated  groups  and  that  of  the  control 
group.  Correlations  between  QuantiGene  and  microarray 
data  were  calculated  using  GraphPad  Prism®  5.04  (Graph- 
Pad  Software,  Inc.,  La  Jolla,  CA). 

The  Pearson  product  moment  (linear)  correlation  coef¬ 
ficient  and  correlation  significance  for  each  transcript  in 
the  microarray  and  Quantigene  bead  platforms  were  cal¬ 
culated  using  base  package  of  R.  Correlation  matrix  was 
calculated  and  plotted  using  corrplot  package  of  R. 

Data  deposition 

All  datasets  used  in  the  study  have  been  deposited  in  Gene 
Expression  Omnibus  of  NCBI  [GEO  accession  #:  GSE45035] 
[http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE45035]. 

Additional  files 


Additional  file  1:  Supplemental  Figures  S1  S20. 

Additional  file  2:  Table  SI.  Pathways  and  process  associated  with 
differentially  regulated  transcripts  across  brain  regions  and  time  points 
(Note:  the  p  value  is  the  gene  level;  but  the  pathways  and  processes 
were  corrected  using  FDR  at  q  <  0.05). 

Additional  file  3:  Table  S2.  Pathways  and  process  associated  with 
differentially  regulated  transcripts  that  passed  p  <  0.01  and  fold 
change  >  2.0  (across  brain  regions  and  time  points)  (Note:  the  p  value 
is  the  gene  level;  but  enrichmenst  of  the  associated  pathways  and 
processes  were  corrected  using  FDR  at  q  <  0.05). 


Abbreviations 

PTSD:  Post  traumatic  stress  disorder;  Agg  E:  Aggressor  exposed  (social 
stressed);  C  Ctrl:  control;  DEG:  Differentially  expressed  gene; 

FIC:  Hippocampus;  AY:  Amygdala;  MPFC:  medial  prefrontal  cortex;  VS:  Ventral 
striatum,  includes  nucleus  accumbens  (NAc);  SE:  Septal  region;  ST:  Corpus 
striatum;  T10  and  T5:  10  days  and  5  days  of  aggressor  exposure  sessions;  Rl, 
RIO  and  R42:  One,  10  and  42  (post  Aggressor  exposure)  days,  respectively; 
GPCR:  G  protein  coupled  receptors;  FDR:  False  discovery  rate. 
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